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Abstract 



A general energy-based theory of autoresonance (self-sustained resonance) 
in low-dimensional nonautonomous systems is presented. The equations that 
together govern the autoresonance solutions and excitations are derived with 
the aid of a variational principle concerning the power functional. These equa- 
tions provide a feedback autoresonance-controlling mechanism. The theory is 
applied to Duffing-like systems to obtain exact analytical expressions for au- 
toresonance excitations and solutions which explain all the phenomenological 
and approximate results arising from a previous (adiabatic) approach to au- 
toresonance phenomena in such systems. The theory is also applied to obtain 
new, general, and exact properties concerning autoresonance phenomena in a 
broad class of dissipative and Hamiltonian systems, including (as a particular 
case) Duffing-like systems. 

PACS number: 05.45.-a Nonlinear dynamics and nonlinear dynamical 
systems 

It has been well known for about half a century that autoresonance (AR) phenomena 
occur when a system continuously adjusts its amplitude so that its instantaneous nonlin- 
ear period matches the driving period, the effect being a growth of the system's energy. 
Autoresonant effects were first observed in particle accelerators [1,2], and have since been 
noted in nonlinear waves [3,4], fluid dynamics [5,6], atomic and molecular physics [7,8], 
plasmas [9-11], nonlinear oscillators [12,13], and planetary dynamics [14-17]. Apparently, 
the first mention of the notion of resonance ( "risonanza" ) was by Galileo [18]. Remarkably, 
this linear-system-based concept has survived up to now: resonance (nonlinear resonance) is 
identified with how well the driving period fits (a rational fraction of) a natural period of the 
underlying conservative system [19]. However, the genuine effect of the frequency (Galilean) 
resonance (FR) (i.e., the secular growth of the oscillation amplitude) can no longer be ob- 
served in a periodically driven nonlinear system. As is well known, the reason is simple: a 
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linear oscillator has a single period which is energy-independent, while nonlinear oscillators 
generally present an infinity of energy- dependent periods. This means that, although an FR 
can still be momentarily induced in a nonlinear system by exciting it with a driving period 
that exactly matches the intrinsic period of the current motion, the subsequent growth of the 
nonlinear oscillations changes the intrinsic period of the motion, which no longer matches 
the excitation period and thus takes the system out of FR. Since linear oscillations represent 
a limiting degenerate (energy) case of the more general nonlinear oscillations, it seems that 
any truly nonlinear generalization of the notion of resonance, in its early (etymological) 
sense of resonare (i.e., awaken an echo of some underlying nonlinear oscillation), should 
be based on energy (or action) considerations. A case has been provided by the notion of 
geometrical resonance [20]. Thus, if one is interested in obtaining a nonlinear equivalent 
of the secular maintained growth intrinsic to the FR, it is clear that the system must not 
be driven by a strictly periodic excitation. In this regard, a previous theoretical approach 
to autoresonance phenomena [3,7-11] provided an early explanation of the mechanism in- 
ducing the growth of the oscillation (without the use of feedback) for particular classes of 
resonantly driven nonlinear systems which stay locked with an adiabatically varying per- 
turbing oscillation (the drive). The adiabatic excitation yields the autoresonant effect by 
automatically adjusting the system's amplitude so that the instantaneous nonlinear period 
matches the driving period. It should be stressed that a fundamental part (hereafter referred 
to as adiabatic autoresonance (AAR) theory, cf. refs. [9-11,21,22]) of the aforementioned 
previous theoretical approach to AR phenomena presents severe limitations of applicability 
and insight: essentially (see ref. [22] for a review), it was developed for nonlinear oscillators 
that reduce to a Duffing oscillator 

x +Uq(x + bx 3 ) = -5 x +e cos (u t + at 2 / 2 ) (1) 

for small amplitudes, where a is the linear sweep rate and 5 > 0. In the context of AAR 
theory, it has been found numerically that AR solutions only occur if (i) the damping 
coefficient 5 is not too large, and (ii) the amplitude of the AR oscillations grows on the 
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average, but also oscillates around the average growth. Also, AAR theory predicts that 
(iii) there exists a threshold for AR, in particular, if the normalized excitation amplitude 
e/to^ 2 exceeds a threshold proportional to a 3 / 4 , the system will follow the excitation to 
high amplitude, while the amplitude will stay very low otherwise, (iv) that the threshold 
sweep rate a th scales as 5 2 , (v) that the AR effect is solely expected for the case with initial 
conditions near some equilibrium of the (unperturbed) nonlinear system, and (vi) that there 
exists a breaking time for AR, t fe . Properties (ii), (iii), (v), (vi) also hold in (vii) the case 
with no dissipation (cf. refs. [9, 10, 21]), but there has as yet been no theoretical explanation 
of that fact. It is worth mentioning that, to the best of the author's knowledge, the case of 
weak dissipation has only been considered in a single previous work (cf. ref. [11]). 

In this Letter a new, general, and energy-based theory for AR phenomena in nonau- 
tonomous systems is presented and applied to the above Duffing oscillators to explain con- 
jointly points (i)-(vii) as well as to deduce new properties concerning AR phenomena in 
generic systems (including Duffing-like systems). The theory arises from the question as to 
whether there exists an upper limit for the growth rate of the system's amplitude when a 
sma//-amplitude force acts on the system. Consider the general family of systems 

x— g(x) — d{x, x) +p(x : x)F(t) : (2) 

where g(x) = — dV(x)/dx [V(x) being an arbitrary time-independent potential], — d(x, x) is 
a general damping force, and p(x, x)F(t) is an as yet undetermined suitable AR-inducing 
force. Clearly, the corresponding equation for the energy is E = x [— d(x, x) + p(x, x)F(t)] = 
P(x, x, t), where E(t) = (1/2) x (t)+V[x(t)] and P(x, x, t) are the energy and power, respec- 
tively. In the spirit of the aforementioned energy-based approach to resonance phenomena, 
the AR solutions are defined by imposing that the energy variation AE = f** P(x,x,t)dt 
is a maximum (with ti,t 2 arbitrary but fixed instants), where the power is considered as a 
functional. This implies a necessary condition (hereafter referred to as the AR condition) 
to be fulfilled by AR solutions and excitations, which is the Euler equation [23] 



From eq. (3), a relationship between x, x, and F can be deduced such that the solutions of 
the system given by eqs. (2) and (3) together provide the AR excitations, -FarW, and the 
AR solutions, XAn(t). It is worth noting that the AR condition (3) represents a feedback 
AR-controlling mechanism, which is absent in the aforementioned previous approach to AR 
phenomena [3,7-11] where an explicit, coordinate-independent, and adiabatic force is used 
from the beginning. In this regard, autoresonant control has been previously discussed in 
the context of vibro-impact systems [24] on the basis of the analysis of nearly sinusoidal self- 
oscillations [25] where the term self-resonance was introduced to indicate "resonance under 
the action of a force generated by the motion of the system itself (cf. ref. [25], p. 166). 
The corresponding AR equations for the multidimensional case can be straightforwardly 
obtained from the same principle and they will be discussed elsewhere [26]. To compare 
the present approach with the previous one [9,11,21,22] (cf. eq. (1)), consider the power 
functional P(x,x,t) —x [—5x +F(t)]. For the particular case of Duffing oscillators, the 
system (2), (3) reduces to 



Note that eq. (4b) gives the AR condition, i.e., the AR excitations and the (corresponding) 
AR solutions have the same instantaneous nonlinear period, at all instants, but without the 
adiabaticity requirement of the AAR theory. Generally, the AR condition (3) means that 
the instantaneous period of the AR solution fits a rational fraction of the instantaneous 
period of the AR excitation. This property generalizes (and contains as a particular case) 
the persistent phase-locking condition of the AAR theory. To obtain AR solutions (and 
hence AR excitations, cf. eq. (4b)) consider the ansatz XAnif) = jf(t) cn \j3g(t) +<f>;m], 
where cn is the Jacobian elliptic function of parameter m, and where the constants (3, m, 
and the functions f(t),g(t) have to be determined for the ansatz to satisfy eq. (4a), while 
7, are arbitrary constants. After substituting this ansatz into eq. (4a), one finds the exact 
general AR solution 



xar { x ar + bx 3 AR ) = 5 x AR , 



(4a) 



Far = 25 x A r ■ 



(4b) 
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x AR (t) =7oe 5t/3 cn [<p (*) ; 1/2] , 
V? (t) = 370^0^ (e 5 */ 3 - 1) /S + ipo, (5) 

with the constraint oo 2 = 25 2 /9 and where <po = <fi + 370^0^/^,70 = 7. Clearly, the exact 
AR excitation corresponding to solution (5) is 

F AR (t) = 2 - l0 5 2 e St ^ cn [<p(t); 1/2] - 2 1 2 5u Vbe 25t / 3 sn [<p (t) ; 1/2] dn [<p (t) ; 1/2] , (6) 

where sn and dn are the Jacobian elliptic functions. Observe that the particular time- 
dependence of the AR solution (5) directly explains the above point (ii) (see fig. 1). 
In comparing the present predictions with those from AAR theory, recall that the lat- 
ter solely exist for the case with x(0) ~ 0, x (0) ~ 0, for b > (point (v)). 
Thus, for this case 70 ~ and hence eq. (6) can be approximated by Far(£) — 
|7 5 2 (l + |5i + ...) cn 70^ (cu t + |cj <5t 2 + ...) ; 1/2 , and, using the Fourier expansion of 
cn [27], one finally obtains 



F AR (t) ~ ^k 70 (J 2 ^1 + ^St + cos Kfj Vb (u t + ^u 5t 2 + 



(7) 



where k ee tt ^2csch(7r/2)/K(l/2) ~ 1, k' = n/(2K(l/2)) ~ 1. Now, one sees that to 
consider the excitation e cos (uj t + at 2 /2) (cf. eq. (1)) as a reliable approximation to Far(£) 
(cf. eq. (7)) implies that the damping coefficient has to be sufficiently small (point (i)) so 
as to have a sufficiently large breaking time, t b ~ 5" 1 (point (vi)). Thus, for t < t b) one 
obtains e th ~ 5 2 ,a th ~ cu 5 (cf. eqs. (1), (7)). When u ~ 5 (recall that tu 2 = 28 2 /9 for 
the exact AR solution (5)), one finds a t h ~ 5 2 (point (iv)), which explains the adiabaticity 
requirement of AAR theory for dissipative systems, e t h/^ 2 ~ cc^/ai 1 ^ 4 = a 3 / 4 (point (iii)), 
and the cosine's argument in eq. (7) can be reliably approximated by the first two terms, 
as in AAR theory (cf. eq. (1)). Figure 1 shows an illustrative comparison between the AR 
responses yielded by AR excitations given by ecos(u)ot + at 2 /2), where in all cases e > e t h, 
and F AR (t) (cf. eq. (6)), respectively, for the cases u ~ 5 (fig. la) and uj 3> 5 (fig. lb). 
Point (vii) is rather striking in view of the very different properties of Hamiltonian and 
dissipative systems, and its explanation is a little more subtle. Firstly, note that current AR 
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theory provides an unsatisfactory result for the limiting Hamiltonian case. For example, eq. 



including the cases of external and parametric (of a potential term) excitations. Clearly, 
the two possible types of corresponding particular solutions, equilibria and those yielded 
by a constant excitation (cf. eqs. (2), (3)), can no longer be AR solutions. Secondly, 
for the above Duffing oscillators one has x AR +u>l (xar + bx AR ) = Far/2 (cf. eq. (4)). 
Therefore, it is natural to assume the ansatz F(t) = X x (t),X > 0, for the case with 
no dissipation, where now the AR rate, A, is a free parameter which controls the initial 
excitation strength. Thus, the corresponding AR solutions are given by eq. (5) while AR 
excitations are given by the expression in eq. (6) multiplied by 1/2, both with A instead 
of S, which explains point (vii) and hence the adiabaticity requirement of AAR theory for 
Hamiltonian systems (recall point (iv)). It is worth mentioning that this valuable result 
holds for the broad family of dissipative systems x +dV(x)/dx = —5 x \x\ + F(t), 
where V(x) is a generic time-independent potential and — 8 x |i| n 1 is a general dissipative 
force (5 > 0,n = 1,2,3,...). The corresponding AR equations (cf. eqs. (2), (3)) are 
x A r +dV(x AR )/dx AR = n5 x AR \x A r\ TI \ F AR = (n + 1)5 x A r \x AR \ n \ and hence one 
obtains x A r ~\~dV (x ar) I dx ar = nFAR/(n + 1). For the limiting Hamiltonian case (5 = 0), it 
is therefore natural to assume the ansatz F(t) = nX x |x| n 1 , A > 0. Thus, AR solutions are 
the same for the dissipative and Hamiltonian cases, while the AR excitations associated with 
the Hamiltonian case are the (corresponding) AR excitations associated with the dissipative 
case multiplied by with A instead of 5 for the Hamiltonian case [28]. In the light 
of the exact AR excitation (cf. eq. (6)), one can readily obtain a reliable approximation 
for arbitrary initial conditions, i.e., not just those near the equilibrium of the unperturbed 
Duffing oscillator: 



where re" = n 2 V2 sech (tt/2) /^ 2 (l/2) ~ 1.61819 ~ (l + y/E) /2 (i.e., the golden ratio). 
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(3) yields r(x) F (t) — for the family given by eq. (2) with d(x, x) = 0,p(x, x) = r(x), i.e., 




(8) 



Thus, for t < t b ~ 5 1 (\ x ) one obtains the general (i.e., valid for any initial condition) 
lst-order adiabatic excitation 

F A ,i(t) = e cos (iu t + at 2 / 2 ) - e ' sin (^o* + at 2 /2) , (9) 

with the above scalings for e th , a t h, and e' th ~ e t hlob 1 / 2 . Figure 2 shows illustrative examples 
for several initial conditions far from x (0) —x (0) = 0. Another fundamental consequence of 
the present approach is the derivation of the scaling laws for the thresholds corresponding to 
higher-order chirps [29]. Indeed, consider the general nth-order adiabatic excitation FA, n (t) = 
e cos [uj (t) t] — e' sin [uj (t) t], uj (t) = uj + Yl^Li a nt n , instead of e cos (uj t + at 2 / 2) in eq. 
(1), where a n is the nth-order sweep rate (cci = a/2). For this general case, the above 
analysis straightforwardly yields the scaling law e th /uj^ 2 ~ iV (n)a^l n+2 ^ for t < t b ~ 
5 _1 (A _1 ), where a n>t h is the threshold nth-order sweep rate and N(n) = [3™ (n + i)!] 3 /( 2n+2 ) 
is a monotonous increasing function. Thus, the 3/4 scaling law is a particular law which 
solely applies to a linear chirp. For the case of a single chirp term (uj (t) = uj + ot n t n ), the 
dependence of the above general scaling law on n indicates that one can expect a similar 
AR effect for ever smaller values of a n as n increases. Computer simulations confirm this 
point: an illustrative example is shown in fig. 3. 

A further question remains to be discussed: We have seen why AAR theory requires 
AR excitations to be adiabatically varying perturbing oscillations, but which are the un- 
derlying adiabatic invariants? To answer this question, note that eq. (4a) (with A instead 
of S for the case with no dissipation) can be derived from a Lagrangian, which one de- 
fines as L = e~ 5t (p 2 /2 — uo 2 x 2 /2 — uo 2 bx 4 /A) , p =x, and whose associated Hamiltonian is 
H = p 2 e St /2 + u 2 (x 2 /2 + bx A /4)e~ 5t . The form of this Hamiltonian suggests the following 
simplifying canonical transformation: X = xe~ St / 2 ,P = pe 5t ^ 2 . It is straightforward to see 
that the generating function of the canonical transformation [30] is F 2 (x,P,t) = xPe~ St l 2 . 
The new Hamiltonian therefore reads: K(X, P, t) = H(x,p, t)-dF 2 /dt = P 2 /2 + lu 2 (X 2 /2 + 
be 5t X 4 /4:) + 5PX/2. In the limiting linear case (6 = 0), one sees that K is conserved, i.e., the 
AR solutions corresponding to the linear system are associated (in terms of the old canoni- 
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cal variables) with the invariant e 5t p 2 /2 + uole~ St x 2 /2 + 5xp/2, while for the nonlinear case 
(b 7^ 0) one obtains (after expanding e 5t ) that the respective AR solutions are associated 
with the adiabatic invariant p 2 /2 + u 2 (x 2 /2 + 6x 4 /4) + 5xp/2 = E + Sxp/2 over the time 
interval ^ t ^ £at, ^a/ ~ (i-e., the same scaling as for the breaking time, tf>, deduced 
above), where E is the energy of the underlying integrable Duffing oscillator. Observe that 
the adiabatic invariant reduces to E provided that 5 (A) is sufficiently small (as required in 
AAR theory) and that the same result is obtained for a general potential V(x) instead of 
Duffing's potential. 

In sum, a general energy-based theory of AR phenomena in low-dimensional nonau- 
tonomous systems has been deduced from a simple variational principle concerning the 
power functional. In particular, the theory explains all the phenomenological and approxi- 
mate results arising from a previous adiabatic approach to AR in Duffing-like systems. For 
this class of systems, the present theory also explains the adiabaticity requirement as well 
as why the same theoretical predictions hold in the cases with and without dissipation, and 
yields the analytical expression for the adiabatic invariants. Additionally, new adiabatic 
approximations to AR excitations are derived concerning two general cases which were not 
considered in the previous adiabatic approach, namely, the case of arbitrary initial condi- 
tions (not just those near equilibria) and the case of arbitrary potential (not just linear) 
chirps, for which new general scaling laws were deduced (including the 3/4 scaling law as 
a particular case). Computer simulations confirmed all the theoretical predictions. In view 
of the generality of the present theory of AR, one can expect it to be quite readily testable 
by experiment (e.g., in the Diocotron system in pure-electron plasmas), and that it will find 
applications in different fields of physics, such as plasmas, fluids, and solar system dynamics. 
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was the origin of the present work, and Professors Gallas, Goldhirsch, Malhotra, Malomed, 
and Sanchez for useful comments. The author acknowledges financial support from Spanish 
MCyT through FIS2004-02475 project. 
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A. Figure Captions 



Figure 1. Autoresonant responses (energy vs time, both variables in arbitrary units) to a 
linearly swept excitation (cf. eq. (1), grey lines) and to an exact AR excitation (cf. eq. (6), 
black lines), for the parameters b = 5,x(0) = lCT 3 ,:r (0) = 0,7 = 10" 3 ,y?o = 0. (a) Case 
with c^o ~ 5, as required for an exact AR excitation, and e = 0.5. (b) Case with u ^> 5, 
i.e., far from the exact AR excitation requirement, and u = 2n. 

Figure 2. Autoresonant responses (energy vs time, both variables in arbitrary units) 
to a general lst-order adiabatic excitation (cf. eq. (9), black lines) and to a harmonic 
and linearly swept excitation (cf. eq. (1), grey lines), for the parameters b = 5, 5 = 
0.4, uj = 0.2, e = 0.5 ~ e t h,a = 0.08 ~ a t h, and initial conditions x(0) = 0.8, x (0) = 0.107 
( 7o = 0.8, e' = 0.9 ~ e' th , thick lines) and x(0) = 0.6, x (0) = 0.08 (70 = 0.6, e' = 0.67 ~ e' th , 
thin lines). 

Figure 3. Autoresonant responses (energy vs time, both variables in arbitrary units) 
to a harmonic excitation with a linear chirp (cu(t) = u + a\t, cf. eq. (1), grey lines) 
and with a quadratic chirp {u (t) = u + a2t 2 , black lines), for the parameters b = 5, 5 = 
0.4, u>o = 0.2, e = 0.5 ~ e t h, ot\ = 0.04 ~ ai,th, «2 = 0.003 ~ a2,th, and the initial conditions 
x(0) = 10~ 3 ,x (0) = (thick lines) and x(0) = 0,x (0) = 1 (thin lines). 
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